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Abstract:  One  of  the  problems  associated  with  the  use  of  the  sodium 
laser  guide  star  as  a  beacon  for  an  adaptive  optics  system  is  the  changes 
in  the  sensitivity  of  the  wavefront  sensor  caused  by  the  extended  nature 
of  the  sodium  layer.  Shack-Hartmann  wavefront  sensors  use  quad  cell 
detectors  to  estimate  the  amount  of  x-  and  y-  gradients  in  each  sub¬ 
aperture.  The  gradients  are  multiplied  by  a  constant  gain  factor  that 
is  equal  in  the  x  and  y  directions.  Due  to  the  non-zero  thickness  of  the 
sodium  layer,  the  wavefront  sensor  sees  an  elongated,  rather  than  round, 
spot,  which  causes  different  gains  in  the  x-  and  y-directions.  This  di¬ 
rectly  impacts  the  reconstructed  phase  and  thus  degrades  the  closed 
loop  AO  performance.  This  paper  presents  an  experimental  demonstra¬ 
tion  of  the  null  space  technique  that  allows  for  real  time  extended  spot 
gain  correction  in  a  closed  loop  AO  operation.  The  implementation  of 
the  technique  in  the  AS  ALT  laboratory  testbed  is  described,  system  per¬ 
formance  degradation  for  incorrect  gains  is  measured,  and  null  space 
corrected  performance  for  different  atmospheres  is  presented. 
Keywords:  adaptive  optics,  wavefront  sensor,  laser  guidestar,  sodium 
beacon,  perspective  elongation 
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1.  Nomenclature 


•  AO  -  Adaptive  Optics 

•  AFRL  -  Air  Force  Research  Laboratory 

•  ASALT  -  Atmospheric  Simulation  and  Adaptive  Optics  Laboratory  Team 

•  ACS  -  Atmospheric  Compensation  Simulator  (software  package) 

•  ATS  -  Atmospheric  Turbulence  Simulator  (hardware  on  the  table) 

•  DM  -  Deformable  Mirror 

•  LGS  -  Laser  Guide  Star 

•  NGAS  -  Sodium  Guidestar  Adaptive  Optics  for  Space  Situational  Awareness 

•  NUGE  -  Non-Uniform  Gain  Experiment 

•  RTR  -  Real  Time  Reconstructor 

•  SH  -  Shack  Hartmann 

•  SM  -  Steering  Mirror 

•  SRI  -  Self  Referencing  Interferometer 

•  WFS  -  Wavefront  Sensor 


2.  Introduction 

The  Air  Force  Research  Laboratory’s  Starfire  Optical  Range  uses  a  sodium  laser  guide 
star  as  a  beacon  for  its  3.5  m  telescope  AO  system.  Sodium  laser  guide  stars  have  been 
used  for  the  purpose  of  beacons  in  adaptive  optics  since  1997  [1].  While  they  provide 
significantly  improved  closed  loop  performance  over  Rayleigh  beacons,  they  also  pose 
numerous  problems.  Some  of  those  problems  are  associated  with  the  sodium  layer,  such 
as  its  varying  thickness,  changing  density  of  the  scatterers  (column  abundance),  and  ex¬ 
tended  nature  of  the  formed  beacon.  This  paper  focuses  on  the  extended  nature  of  the 
beacon  and  its  effects  on  the  LGS  WFS  response,  which  leads  to  the  degradation  of  the 
reconstructed  phase,  and  consequently  decreased  Strehl  ratio. 

The  mesospheric  sodium  layer  is  located  at  about  90  km  above  the  sea  level,  and  has  a 
thickness  that  varies  from  8  to  16  km  [2].  Due  to  this  considerable  thickness,  a  LGS 
created  in  this  layer  by  a  laser  beam  tuned  to  the  sodium  Di  line  at  589  nm,  appears  as  an 
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extended  source  (see  Figure  1).  The  extended  nature  of  the  source  presents  two  potential 
problems  for  WFS.  First,  the  spots  will  appear  extended  on  the  WFS  camera.  Second, 
the  spots  will  look  different  from  sub-aperture  to  sub-aperture,  and  thus  the  WFS  gains 
in  each  sub-aperture  will  differ  from  one  another. 

Numerous  models  have  been  proposed  by  the  astronomical  community  that  estimate  the 
effects  of  both  the  extended  nature  of  the  beacon  and  the  per-sub-apeilure  gain  varia¬ 
tions  [3],  [4],  [5).  Most  common  techniques  include  the  per-sub-aperture  specific  gains 
applications  [6],  noise-optimal  matched  filter  spot  position  estimation  algorithm  [7],  [8], 
dynamic  refocusing  of  the  telescope  [9],  and  dynamically  steered  SH  sub-apertures  [10]. 

In  this  paper  we  present  a  perspective  elongation  solution,  based  on  the  slope  discrepancy 
space  of  the  reconstructor  matrix.  Perspective  elongation  is  assumed  to  vary  only  in  the 
x-  and  y-  directions  of  the  primary  aperture,  and  not  on  the  sub-aperture  basis.  The  null 
space  of  the  reconstructor  is  used  in  order  to  separate  extended  beacon  effects  from  the 
atmospheric  turbulence,  reference  projection  into  the  null  space  is  compared  with  the 
frame-by-frame  gradients,  and  correction  factor  is  calculated  and  applied  to  gradients  in 
real  time. 


3.  Effects  of  the  Extended  Source 

Figure  1  shows  the  geometry  of  the  sodium  layer  with  relation  to  the  primary  imaging 
telescope,  and  the  sodium  beacon  launching  telescope.  Distance  between  the  launching 
and  imaging  telescopes,  thickness  of  the  sodium  layer,  and  the  zenith  angle  are  the  pri¬ 
mary  variables  responsible  for  the  extended  nature  of  the  sodium  LGS. 
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Fig.  1 .  Geometry t  of  the  sodium  layer,  Na  guide  star,  and  the  telescope. 
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In  astronomical  AO  systems,  Shack-Hartmann  WFSs  are  most  commonly  used  for  sens¬ 
ing  the  wavefront  of  LGS.  The  sensor  consists  of  a  lenslet  array,  where  each  lens  (each 
sub-aperture)  is  followed  by  a  quad  cell  detector  (2x2  array  of  active  pixels),  or  higher 
resolution  detector,  which  is  typically  binned  down  to  a  2  x  2  quad  cell.  A  typical  quad 
cell  with  a  focused  spot  on  it  is  shown  in  Figure  2.  The  gradients  for  such  a  sub-aperture 
are  calculated  using  formulas  I  and  2,  where  a,  b,  c,  and  d  are  light  intensities  measured 
by  corresponding  pixels,  and  subsequently  multiplied  by  a  calibration  factor,  which  takes 
into  account  local  measured  biases. 


x grads  = 
y grads  — 


(, b  +  d)-(a  +  c ) 
a+b+c+d 
(a  +  b)  -  ( c  +  d ) 
a+b+c+d 


(1) 

(2) 


Sensitivity  of  the  WFS  for  a  particular  spot  size  can  be  determined  by  scanning  the 
spot  across  the  quad  cell,  and  plotting  calculated  gradients  versus  the  spot  centroid. 
The  resulting  plot  is  called  a  transfer  curve  (or.  sometimes,  the  S-  curve),  where  the 
slope  in  the  middle  is  the  sensitivity  of  the  WFS,  and  the  inverse  of  the  slope  is  the 
gain.  For  a  round  spot,  x-  and  y-  gains  are  equal  (see  Figure  2).  Larger  spots  pro¬ 
duce  less  steep  curves,  the  sensitivity  goes  down,  and  the  gains  of  the  WFS  go  up, 
as  is  illustrated  in  Figure  3.  When  the  object  being  imaged  is  extended,  the  quad 
cell  has  different  x-  and  y-  responses,  which  results  in  the  transfer  curves,  and  sub¬ 
sequently  sensitivities  and  gains  being  different  in  the  x-  and  y-  directions.  Figure  4 
shows  a  simulated  oblong  spot  along  with  corresponding  x-  and  y-  transfer  curves. 


Transfer  curvet  for  small  spot 


spot  centroid 


Fig.  2.  Small  round  spot  at  the  center  of  the  quad  cell  produces  identical  x- 
and  y-  transfer  cun  es,  and  thus  equal  x-  and  v-  gains. 
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Transfer  curves  for  small  spot 


Fig.  3.  Large  round  spot  produces  equal  transfer  curves  in  x-  and  y-  di¬ 
rection,  but  the  transfer  curves  are  different  from  the  small  spot  curves. 
In  particular,  the  slope  in  the  middle  is  less  for  larger  spot,  resulting  in 
greater  required  WFS  gains. 


Transfer  curves  In  x-  and  y-  for  oblong  spot 


a 

b 

c 

d 

Fig.  4.  Oblong  spot  produces  different  quad  cell  responses  in  x-  and  y- 
directions.  This  results  in  different  transfer  curves  with  different  slopes. 
The  x-  and  y-  WFS  gains  are  different. 


Figure  5  shows  ACS  software  simulation  data  for  a  high  resolution  SH  camera  for 
sodium  beacon  LGS  with  mild  turbulence.  The  apparent  extension  of  the  spot  varies  from 
sub-aperture  to  sub-aperture,  though  only  one  sub-aperture  is  shown.  Produced  transfer 
curves  for  x-  and  y-  directions  are  shown  on  the  right  in  Figure  5.  It  is  important  to  point 
out  that,  even  though  the  curves  look  almost  identical,  the  change  in  sensitivity,  and  thus 
the  change  in  gains  between  the  x-  and  y-  directions  is  significant.  Data  calculated  di¬ 
rectly  from  Figure  5  is  shown  in  Table  1,  where  the  change  in  going  from  x-  to  y-  gains  is 
14.54%.  This  is  significant  and  will  noticeably  degrade  the  Strchl  ratio  of  the  closed  loop 
performance.  The  problem  that  this  paper  addresses  is  dynamically  changing  gains  of  the 
WFS,  which  are  caused  by  the  extended  nature  of  sodium  beacon,  and  are  uncorrected  in 
normal  closed  loop  AO  runs. 
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Transfer  curves  for  ACS  NaB  simulation 


Fig.  5.  ACS  simulated  sodium  beacon  spot,  one  sub-aperture,  with  mild 
turbulence. 


Sodium  Beacon,  1  sub-aperture 

X  sensitivity 

0.6781 

Y  sensitivity 

0.7935 

X  gain 

1.4748 

Y  gain 

1.2603 

Table  1.  Sensitivities  and  gains  in  x-  and  y-  for  NaB  ACS  simulated  spots, 
one  sub-aperture. 

4.  Theoretical  Solution  -  The  Null  Space  of  the  Reconstructor 

The  solution  to  perspective  elongation  presented  in  this  paper  identifies  system  degra¬ 
dations  due  to  the  atmosphere  and  elongated  sources  separately.  AFRL  developed  the 
Nullspace  algorithm  for  gain  correction,  which  uses  the  null  space  of  the  reconstructor  to 
look  for  gradient  bleed  through.  The  null  space  (slope  discrepancy  space)  of  the  recon¬ 
struction  matrix  contains  very  little  reconstructible  atmospheric  energy  (almost  none), 
however  it  docs  see  the  gradient  components.  By  projecting  the  reference  gradient  com¬ 
ponents  into  the  null  space,  and  comparing  them  with  frame-by-frame  gradient  projec¬ 
tions,  we  calculate  the  so-called  gradient  "bleed  through”,  which  shows  the  effects  of  the 
extended  beacon  on  the  frame-by-frame  gradients. 


4. 1.  Null  Space  of  the  Reconstructor 

To  explain  the  concept  of  null  space  in  more  detail,  we  start  with  the  G  matrix,  which  is 
used  to  map  the  phase  into  gradients.  Figure  6  (left)  shows  the  G  matrix  decomposed  into 
three  separate  matrices  by  singular  value  decomposition  method:  G  =  U  ■  S  -V 7  .  The  V 
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matrix  consists  of  the  orthogonal  unit  vectors  that  span  the  space  of  active  mirror  actua¬ 
tors.  The  U  matrix  consists  of  orthogonal  unit  vectors  that  span  the  space  of  all  possible 
gradients.  Note  that  the  number  of  active  gradients  is  greater  than  the  number  of  actua¬ 
tors  by  about  a  factor  of  two,  which  means  the  system  is  overdetermined.  The  S  matrix 
takes  the  amount  of  one  mode  in  V  and  scales  it  to  the  corresponding  element  of  U.  By 
taking  the  pseudo-inverse  of  G,  the  RTR  matrix  is  obtained,  which  takes  in  the  gradients 
and  reconstructs  the  phase.  The  RTR  matrix  can  be  expressed  similarly  in  terms  of  U, 
S,  and  V  matrices:  RT R  =  V  ■  Smv  ■  UT  -  see  Figure  6  (right).  The  rows  of  the  U1  matrix 
span  the  space  of  all  possible  gradients  (it  creates  the  gradient  "modes”).  By  multiplying 
UT  by  the  vector  of  frame-by-frame  gradients,  we  get  the  amount  of  gradients  in  each  of 
the  possible  gradient  modes.  S,„v  matrix  converts  possible  gradients  to  mirror  actuators. 
However,  only  the  first  M  x  2N  gradient  modes  (where  M  is  the  number  of  active  actua¬ 
tors,  2N  is  the  number  of  possible  gradients)  are  converted  into  the  reconstructed  phase, 
which  corresponds  to  possible  reconstructible  modes  of  the  reconstructor.  The  remaining 
[2/V  -  M]x2N  modes  are  multiplied  by  zero,  they  are  not  affected  by  the  atmospheric 
energy  and  are  not  manifested  in  the  resulting  reconstructed  phase.  Those  modes  span  the 
null  space.  Note  that  the  null  space  model  was  originally  developed  for  NGAS  program, 
and  the  dimensions  of  the  G,  U,  S,  V  and  RTR  matrices  shown  in  Figure  6  correspond  to 
the  system  parameters  of  NGAS.  Very  little  reconstructible  atmospheric  energy  gets  into 
the  null  space,  with  the  exception  of  high  frequency  modes,  which  cannot  be  corrected 
by  the  DM.  By  looking  at  the  power  spectral  density  of  each  mode  for  reconstructible 
modes  and  for  null  space  modes,  it  is  apparent  that  the  null  space  does  not  contain  much 
of  the  reconstructible  atmospheric  energy.  Figure  7  illustrates  the  concept.  For  the  recon¬ 
structible  modes  (left  plot)  going  from  open  loop  to  closed  loop  (from  blue  to  red  curves), 
the  energy  at  lower  frequencies  decreases  significantly.  The  energy  at  higher  frequencies 
is  unchanged  due  to  the  inability  of  the  closed  loop  AO  to  compensate  for  higher  fre¬ 
quency  turbulence.  By  contrast,  the  power  in  the  null  space  modes  (right  plot)  remains 
near  identical  in  open  (blue  curve)  and  closed  (red  curve)  loop.  This  signifies  the  fact  that 
the  DM  does  not  affect  the  energy  in  null  space,  and  thus  it  can  be  used  to  sense  changes 
in  the  system  unaffected  by  the  atmosphere. 
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Fig.  6.  G  matrix  maps  the  phase  into  the  gradients.  Shown  is  singular  value 
decomposition  of  the  G  matrix.  RTR  matrix  is  obtained  by  pseudo-inverting 
the  G  matrix.  RTR  converts  gradients  into  phase. 
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Fig.  7.  Power  in  the  reconst rucuible  modes  (left)  and  nail  space  inodes 
(right)  versus  frequency  for  open  loop  (blue)  and  closed  loop  (red). 


4.2.  Calculations  behind  the  algorithm 

Applying  the  Nullspace  algorithm  to  closed  loop  AO  data  includes  two  parts  -  the  offline 
part,  which  is  done  once,  and  the  online  part,  which  is  calculated  on  each  frame  of  the 
new  gradients  data.  For  the  new  reconstructor  and/or  new  reference,  the  offline  part  needs 
to  be  calculated  again.  Suppose  the  number  of  active  X-  gradients,  equal  to  the  number 
of  active  Y-  gradients,  is  N;  the  number  of  active  actuators  is  M. 

Offline: 


•  Singular  Value  Decomposition  is  used  to  calculate  the  null  space  of  the  reconstruc¬ 
tor,  this  matrix  is  called  Z  ([2 N  -  M]  x  2 N).  Z  maps  the  gradients  into  the  null 
space. 

•  Zx  and  Zy  vectors  are  formed,  where  Zx  is  the  Z  matrix  applied  to  X-  gradient  com¬ 
ponents  of  the  reference  (reference  with  Y-  components  zeroed),  and  Zy  is  a  vector 
formed  by  applying  the  Z  matrix  to  Y-  components  of  the  reference  (reference  with 
X-  components  zeroed).  Both  Zx  and  Zy  have  dimensions  of  [2/V  —  M\  x  1. 

•  Zx  and  Zy  are  combined  into  one  matrix  [2/V  -  M]  x  2  which  is  pseudo-inverted  to 
get  ModeslSensitivity  matrix,  with  dimensions  2  x  [2/V  -  M ].  This  matrix  maps  the 
null  space  spectrum  of  the  averaged  gradients  into  sensitivity  values  in  X-  and  Y-. 
It  shows  how  much  of  the  X-  and  Y-  reference  is  present  in  the  null  space. 

•  Gradients2Sensitivity  operator  is  found  by  multiplying  Modes2Sensitivity  matrix 
(2  x  [2 /V  —  M\)  by  Z  matrix  ([2/V  —  A7]  x  2N).  The  resulting  matrix  has  dimensions 
of  2  x  2/V,  it  maps  gradients  directly  into  the  sensitivities.  In  other  words,  it  com- 
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bines  the  calculation  of  the  null  space  components  of  the  gradients  with  that  of 
extracting  the  sensitivities  from  those  null  space  components. 

Online: 

•  Gradients2Sensitivity  operator  is  used  on  frame-by-frame  gradients  to  calculate  the 
new  sensitivities. 

•  Gain  is  calculated  as  the  inverse  of  the  sensitivity:  Gain  =  1  /Sensitivity. 

•  The  gains  are  updated  according  to  the  control  law: 

Filtered -gain(i  4-  1 )  =  A  x  leak-gain(i)  +Bx  servo  ^>ain(i).  (3) 


5.  Experimental  Setup 

5.7.  ASALT  Lab  Optical  Bench  Layout 

The  ASALT  Laboratory  provides  a  complete  bench-top  AO  system  with  a  capability 
to  simulate  a  number  of  atmospheric  conditions.  The  layout  of  the  optical  bench  used 
to  test  the  Nullspace  algorithm  is  shown  in  Figure  8.  Note  that  the  system,  along  with 
all  the  atmospheric  parameters  simulated,  is  scaled  to  mimic  a  1 .5  m  telescope  system, 
but  could  also  be  scaled  to  other  telescope  sizes.  A  1550  nm  laser  is  used  to  simulate 
the  point  source  (see  Section  5.2  for  simulation  of  the  extended  source)  that  takes  the 
light  through  the  atmospheric  turbulence  simulator,  which  includes  two  relay  telescopes 
and  two  phase  wheels  etched  with  the  Kolmogorov  structure  function.  By  exchanging 
the  field  lenses  and  positions  of  relay  telescopes,  and  by  controlling  the  speed  of  phase 
wheel  rotations,  atmospheric  parameters  such  as  Fried’s  coherence  length,  Rytov  number, 
and  Greenwood  frequency  can  be  controlled.  Next,  the  aberrated  light  beam  reflects  off 
the  steering  mirror,  which  removes  global  tip/tilt.  The  light  is  then  reflected  onto  off- 
axis  parabolic  mirrors,  bringing  the  beam  to  a  collimated  space,  where  the  DM  applies 
high  order  correction  to  the  beam.  A  Xinetics  941  actuator  DM  with  31x31  actuators 
was  used  for  this  experiment.  After  the  DM,  the  light  is  sent  to  beamsplitters,  which 
separate  the  beam  into  different  paths  -  a  SH  WFS  path,  a  Scoring  Camera  path,  and 
a  Self-Referencing  Interferometer  path.  The  wavefront  sensor  used  for  the  Nullspace 
experiment  was  a  24x24  sub-aperture  Shack-Hartmann  WFS,  with  2x2  active  pixel  quad 
cells.  A  128x128  pixel  scoring  camera  was  used  to  monitor  system  performance. 
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ATS 


1550  nm  source 


SM 


Fig.  8.  ASALT  Lab  SRI-1  table  layout. 


5.2.  Modeling  the  Extended  Source 

While  the  technique  described  in  this  paper  focuses  on  correcting  the  gains  of  the  WFS 
due  to  the  extended  nature  of  the  source,  at  the  time  of  this  experiment,  the  ASALT 
Lab  did  not  have  the  extended  optical  source  to  accurately  mimic  the  sodium  beacon.  In 
order  to  simulate  the  extended  source,  the  gains  of  the  WFS  were  intentionally  degraded 
through  data  post-processing.  Gain  error  was  introduced  by  multiplying  closed  loop  no 
turbulence  (vacuum  propagation)  gradient  gains  by  a  constant  factor  prior  to  closing  the 
loop.  Two  incorrect  gain  factors  were  used:  kx  for  x  gradients,  and  ky  for  y  gradients. 
The  gradients  were  multiplied  as  follows: 

x.  gradients  =  kx  x  x. gradients  (4) 

y.  gradients  =  ky  x  y. gradients  (5) 

Figure  9  shows  the  perfect  spot  for  closed  loop  AO  run  on  the  left,  and  the  spot  created 
by  setting  kx=0.4  and  ky=0.2  on  the  right.  By  varying  the  magnitude  of  kx  and  ky,  the 
elongation  of  the  source  can  be  varied  in  two  different  directions.  For  most  of  the  tests 
presented  in  this  paper,  gain  factors  were  set  at  kx=0.4,  ky=0.2.  Such  values  result  in  poor 
closed  loop  performance  when  the  extended  nature  of  the  source  is  uncorrected.  While 
for  the  real  sodium  beacon  the  performance  degradation  due  to  the  given  thickness  of  the 
sodium  layer  is  nowhere  near  as  dramatic,  the  purpose  of  this  paper  is  to  demonstrate  that 
the  Nullspace  algorithm  works  well  for  the  worst  possible  scenario. 
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Fig.  9.  Left:  Scoring  camera  image  for  vacuum  case,  kx=I.O,  ky-1.0. 
Right:  Scoring  camera  image  for  vacuum  case  for  kx-0.4,  ky-0.2  -  in¬ 
tentionally  set  wrong  gains  create  the  elongation  of  the  source. 


5.3.  Control  System 

The  ASALT  Lab  control  system  uses  an  interactive  environment  that  combines  hard¬ 
ware  and  software  elements,  which  are  run  from  different  computers,  as  components  on 
the  central  layout.  Figure  10  shows  a  snap  shot  of  the  control  layout,  which  captures  one 
frame  (one  iteration)  of  data.  Individual  components  (blocks  in  Figure  10)  are  either  hard¬ 
ware  elements  (DM,  cameras,  tracking  controllers,  ATS),  or  software  elements,  called 
’’processing  nodes”  that  perform  specific  calculations/data  manipulation.  Each  compo¬ 
nent  has  a  certain  number  of  inputs,  outputs,  and  so  called  configuration  items,  allowing 
various  parameters  to  be  changed  dynamically  during  data  runs.  Each  experiment  has  its 
own  layout  designed  for  it,  with  particular  components  built  into  it.  The  flow  of  data  col¬ 
lection  is  established  by  connecting  all  the  required  components  in  a  specified  sequence. 


Fig.  10.  Control  Console,  layout  that  was  used  for  Nullspace  algorithm. 
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5.4.  Parameter  Space 


Table  2  shows  several  sets  of  atmospheric  parameters  that  were  used  to  test  the  Nullspace 
algorithm.  Note  that  for  three  test  runs,  a  combination  of  ro  and  Rytov  were  varied,  rang¬ 
ing  from  severe  to  mild  turbulence.  For  three  other  runs,  ro  and  Rytov  were  held  con¬ 
stant,  but  Greenwood  frequency  was  varied.  For  separate  three  runs,  atmospheric  condi¬ 
tions  were  kept  constant,  but  extended  beacon  gains  kx  and  ky  were  changed  in  order  to 
demonstrate  that  the  algorithm  works  well  with  sources  of  different  perspective  elonga¬ 
tion. 


Data  Set 

Turbulence  Set 

ro(cm) 

Rytov 

fc(Hz) 

kx 

ky 

1 

TR2 

7.8 

0.25 

33.8 

0.4 

0.2 

2 

TR1 

10.3 

0.06 

33.8 

0.4 

0.2 

3 

TR3 

15.32 

0.1 

33.8 

0.4 

0.2 

4 

TR4 

15.32 

0.1 

52.9 

0.4 

0.2 

5 

TR5 

15.32 

0.1 

82.6 

0.4 

0.2 

6 

TR3 

15.32 

0.1 

33.8 

0.4 

0.2 

7 

TR3 

15.32 

0.1 

33.8 

1.5 

0.5 

8 

TR3 

15.32 

0.1 

33.8 

1.8 

1.0 

Table  2.  Parameter  space  for  testing  Nullspace  algorithm  in  AS  ALT  Lab 


6.  Results  and  Analysis 

Strehl  ratio  was  one  of  the  primary  metrics  used  to  analyze  the  performance  of  the  sys¬ 
tem  under  different  atmospheric  conditions.  For  this  experiment,  instantaneous  (frame  by 
frame)  Strehl  was  calculated  as  a  3x3  bucket  power  around  the  hottest  pixel,  normalized 
by  3x3  bucket  power  for  a  vacuum  (closed  loop  no  turbulence  case).  The  left  plot  in  Fig¬ 
ure  1 1  shows  the  instantaneous  Strehl  ratio,  with  the  blue  curve  showing  the  open  loop 
performance,  and  the  pink  curve  showing  the  standard  closed  loop  performance  under 
the  same  atmosphere  without  any  incorrect  gains.  The  black  curve  was  produced  with 
kx=0.4,  and  ky=0.2,  without  correcting  for  those  gain  errors,  which  resulted  in  perfor¬ 
mance  degradation  all  the  way  down  to  open  loop  level.  The  green  curve  shows  the  data 
corrected  by  the  Nullspace  algorithm,  which  started  at  frame  50.  It  is  evident  that  it  took 
about  25  frames  for  the  Strehl  to  start  going  up  significantly.  At  about  frame  150  the 
Nullspace-corrected  data  starts  fluctuating  around  the  regular  closed  loop  AO  data.  Note 
that  a  temporal  control  law  under  the  Nullspace  function  sets  up  how  quickly  the  data  is 
being  updated.  The  right  plot  in  Figure  1 1  shows  the  green  Nullspace-corrected  Strehl 
data  on  top,  and  frame-by-frame  kx  and  ky  gradient  gains  on  the  bottom  (blue  and  red). 
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The  gains  started  at  kx=0.4,  ky=0.2,  as  was  intentionally  set  up  to  simulate  the  extended 
source,  remained  constant  until  frame  50,  and  started  climbing  up  until  they  were  nearly 
1  (correct  value)  for  frames  150  on.  Figure  12  shows  the  frame-by-frame  normalized 


Fig.  1 1.  Instantaneous  Strelil  for  r<)  =  8.87  cm,  Ry  =  0.01,  fc,  =  24.0 
Hz.  Nullspace  corrected  data  (blue)  vs  uncorrected  data  (black).  Right: 
Nullspace  corrected  instantaneous  Streld  (top),  frame-by-frame  gradient 
gain  factors  Lx  and  ky  (bottom). 


power  data  for  the  same  atmospheric  conditions  as  were  used  in  Figure  1 1,  with  the  dif¬ 
ference  that  the  power  in  a  3x3  bucket  was  normalized  by  that  of  the  closed  loop  correct 
gains  case.  By  doing  so.  performance  improvements  due  to  the  Nullspace  algorithm  were 
isolated  from  other  changes  that  might  have  occurred  in  the  system.  Note  that  the  pink 
curve  is  now  always  at  I ,  since  it  is  normalized  by  itself  at  each  frame.  The  blue  and 
green  curves  are  two  identical  Nullspace  runs,  which  start  out  at  the  bottom,  and  end  up 
fluctuating  around  the  best  closed  loop  AO  performance.  There  was  virtually  no  differ¬ 
ence  between  the  two  Nullspace  runs;  they  were  taken  in  order  to  verify  the  repeatability 
of  the  algorithm. 


i  *  •«  *  H3  at  B 


Fig.  1 2.  Data  normalized  by  the  correct  gains  closed  loop  AO.  The  two 
Nullspace  runs  are  identical,  showing  the  repeatability  of  the  algorithm 
performance. 
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Figure  13  shows  several  more  Strehl  plots  that  illustrate  how  well  the  algorithm  worked 
for  the  different  atmospheres  presented  in  Table  2.  Note  that  for  the  data  shown  in  the 
upper  left  corner  in  Figure  13,  we  used  a  different  set  of  field  lenses  compared  to  all  other 
data  sets.  This  could  explain  the  fact  that  this  data  set  has  more  stable  plots  and  higher 
mean  Strehl.  For  data  shown  in  the  lower  left  corner  of  Figure  13,  Strehl  ratio  is  varying 
highly  for  both  correct  and  incorrect  gains.  This  is  due  to  the  fact  that  our  system  was  not 
fast  enough  to  correct  for  high  Greenwood  frequencies.  The  important  point  is  that  the 
uncorrected  extended  gains  plot  is  way  below  the  correct  gains  one  (the  black  curve  is 
way  below  the  pink  curve),  and  the  Nullspace  corrected  plot  lines  up  very  well  with  the 
regular  closed  loop  AO  one.  For  lower  right  plot  in  Figure  13,  the  Greenwood  frequency 
was  set  even  higher,  at  82.6  Hz  (based  on  1  kHz  system  speed).,  however  the  control 
filter  gain  was  adjusted  so  that  the  system  could  keep  up  with  the  turbulence  better.  This 
explains  better  closed  loop  AO  performance  (pink  curve)  compared  to  lower  left  plot 
data,  where  fc  was  52.9  Hz,  but  the  control  gain  was  unadjusted.  The  important  point 
is  that  the  Nullspace  performance  is  consistently  good  for  both  data  sets,  bringing  the 
Strehl  back  to  the  regular  AO  performance  level  in  each  case. 


Co  =  7.8cm:  Ry  =  0.25,  fc  =  33.8//: 


SC  XL 


r0  =  15 .32 cm;  Ry  0. 1 .  fc  =  33.8//: 


ru  =  15.32cm;  Ry  =  0.1.  fc  -  52.9 Hz 


r0  =  15.32 cm;  Ry  =  0. 1 .  fc  =  82.6 Hz 


Fig.  13.  Strehl  Ratio  plots  for  several  atmospheric  sets.  Note  the  Nullspace 
corrected  data  (green)  climbs  up  to  the  closed  loop  AO  case. 
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7.  Conclusion 


We  demonstrated  that  real  time  WFS  gain  correction  technique  based  on  null  space  of 
the  reconstructor  works  well  to  improve  the  closed  loop  AO  performance  in  various  tur¬ 
bulence  regimes.  The  technique  is  relatively  insensitive  to  atmospheric  conditions:  for 
low  to  moderate  turbulence,  the  Nullspace  algorithm  brings  the  system  performance  all 
the  way  up  to  the  best  possible  closed  loop  AO  performance;  for  more  severe  turbu¬ 
lence  cases,  the  Nullspace  technique  does  not  go  quite  all  the  way  to  perfect  closed  loop 
AO  level,  but  exhibits  very  significant  increase  in  Strehl  compared  to  data  uncorrected 
for  extended  beacon  gain  changes.  The  algorithm  produced  stable,  repeatable  data  with 
consistently  improved  Strehl.  In  the  future  we  plan  to  incorporate  a  new  optical  source 
in  the  ASALT  laboratory  that  would  simulate  a  real  extended  beacon  under  different 
zenith  angles.  This  will  allow  changing  gains  in  the  optical  domain,  as  opposed  to  doing 
this  in  software  as  was  done  for  the  experiment  described  in  this  paper.  In  addition,  the 
Nullspace  technique  needs  to  be  combined  with  the  algorithm  that  will  vary  gains  from 
sub-aperture  to  sub-aperture  based  on  the  zenith  angle  and  elevation.  This  will  allow  for 
yet  a  better  method  of  compensation  for  extended  nature  of  the  sodium  beacon. 
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